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I discuss algorithms for the evaluation of Feynman integrals. These algorithms 
arc based on Hopf algebras and evaluate the Feynman integral to (multiple) 
polylogarithms. 



1 Introduction 

Multiple polylogarithms arc an object of interest not only for mathematicians, 
but also for physicists in the domain of particle physics. Here, I will discuss 
how they occur in the calculation of Feynman loop integrals. The evaluation 
of these integrals is an essential part to obtain precise theoretical predictions 
on quantities, which can be observed in experiment. These predictions have a 
direct impact on searches for signals of "new physics" . 

Due to the complexity of these calculations, computer algebra plays an es- 
sential part. This in turn requires that methods developed to evaluate Feyn- 
man loop integrals are suitable for an implementation on a computer. The 
focus for the practitioner shifts therefore from the calculation of a particular 
integral to the development of algorithms for a class of integrals. This shift is 
accompanied by a movement from concrete analytical methods towards ab- 
stract algebraic algorithms. In this article I review some algebraic techniques 
to solve Feynman loop integrals. 

In the next section I show briefly why loop corrections are needed for to- 
day's experiments and how Feynman loop integrals arise in a pratical context. 
Scct.[3]mtroduccs the algebraic tools: A particular form of nested sums, which 
form a Hopf algebra and which admits as additional structure a conjugation 
and a convolution product. In Sect. 0] I show how these tools are used for 
the solution of a simple one-loop integral. Special cases of the nested sums 
are multiple polylogarithms, which are discussed in Sect. These multiple 
polylogarithms admit a second Hopf algebra structure. In Sect. El I discuss 
how the antipode of the Hopf algebra can be used to simplify expressions. 
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2 Phenomenology 

Phenomenology is the part of theoretical particle physics, which is most closely 
related to experiments. The standard experiment in particle physics acceler- 
ates two particles, brings them into collision and observes the outcome. Ex- 
amples for such experiments are Tevatron at Fermilab in Chicago, HERA at 
DESY in Hamburg or LEP and the forthcoming LHC at CERN in Geneva. 
Quite often a bunch of particles moving in one direction is observed. This 
is called a (hadronic) jet and the direction of movement follows closely the 
direction of the initial QCD partons (e.g. quarks and gluons), which where 
created immediately after the collision. The observed events can be classified 
according to their experimental signature, like the number of jets seen within 
one event. Interesting questions related to these events are for example: How 
often do we observe events with two or three jets ? What is the angular dis- 
tribution of these multi-jet events ? What is the value for specific observables 
like the thrust, defined by 

E ll* • n l 

T = max n -^= ri — — , (1) 

i 

which maximizes the total longitudinal momentum of all final state particles 
Pi along a unit vector n ? There are many more interesting observables, and 
theoreticians in phenomenology try to provide predictions for those. Obvi- 
ously, it is desirable not to start a new calculation for each observable, but to 
have a generic program, which provides predictions for a wide class of observ- 
ables. This forces us to work with fully differential quantities. This leaves in 
the main part of the calculation many kincmatical invariants, which cannot 
be integrated out. As a general rule, the higher the number of jets is in an 
event, the more complicated it is to obtain a theoretical prediction, since the 
number of independent kinematical invariants increases. 

By comparing the measurements with theoretical predictions one can de- 
duce information on the original scattering process. By counting the number 
of events with a particular signature one may discover new effects and new 
particles. However, it quite often occurs that hypothetical new particles lead 
to the same signature in the detector as well known processes within the Stan- 
dard Model of particle physics. To distinguish if a small measured excess in 
one observable is due to "new physics" or only due to known physics requires 
precise theoretical prediction from theory In the case where the expected 
events due to "new physics" are only a fraction of the events from Standard 
Model physics, it is most important to have a small theoretical uncertainty for 
the Standard Model background processes. To give a simple example, if the 
number of events due to "new physics" is about 1% of the number of events 
for background processes within the Standard Model, then wc need a theo- 
retical prediction of the background with a precision better than 1%. At the 
energies where the experiments are done, all coupling constants are small and 
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perturbation theory in the coupling constants is the standard procedure to 
obtain theoretical predictions. Therefore to reduce the theoretical uncertainty 
on a prediction requires us to calculate higher orders in perturbation theory. 
In phenomenology we are now moving towards fully differential next-to-next- 
to-leading calculations, e.g. predictions which include the first three orders in 
the perturbative expansion. An essential ingredient of higher order contribu- 




tions are loop amplitudes. Fig. ^ shows a Feynman diagram contributing to 
the one-loop corrections for the process e + e~ — > qgq. From the Feynman rules 
one obtains for this diagram (ignoring coupling and colour prefactors): 

- v(p A )Yv{ Pb ) j±- J * J^ ^ufatffa^lv^l^Yufa)- (2) 

Here, pu =Pi+P2, P123 = Pi +P2 +P3, k 2 = h - P i 2 , k 3 = k 2 -p3- Further 
£(P2) — 1t£ t (p2), where £ T (p2) is the polarization vector of the outgoing 
gluon. All external momenta are assumed to be massless: pf = for i = 1..5. 
Dimensional regularization is used to regulate both ultraviolet and infrared 
divergences. In (jSJ) the loop integral to be calculated reads 

f d^h fc£fcf 

J klklkf {6) 

This loop integral contains the loop momentum k\ in the numerator. The 
further steps to evaluate this integral are now: 

Step 1 Eliminate powers of the loop momentum in the numerator. 
Step 2 Convert the integral into an infinite sum. 
Step 3 Expand the sum into a Laurent series in e. 

The first two steps are rather easy to perform and convert the original integral 
to a more convenient form. The essential part is step 3. It should be noted 
that the integral in © is rather simple and can be evaluated by other means. 
However, I would like to discuss methods which generalize to higher loops and 
this particular integral should be viewed as a pedagogical example. 

To eliminate powers of the loop momentum in the numerator one can trade 
the loop momentum in the numerator for scalar integrals (e.g. numerator = 1) 
with higher powers of the propagators and shifted dimensions E] : 
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2\"3 ■ 



(4) 



This algorithm introduces temporarily Schwinger parameters together with 
raising and lowering operators and expresses one integral of type © in terms 
of several integrals of type . 

In the second step an integral of type is now converted into an infinite 
sum. Introducing Feynman parameters, performing the momentum integra- 
tion and then the integration over the Feynman parameters one obtains 



-)2m-2e 



fcl 



n m-s (-fc 2 )^l (-fcf)^ (-fc|)^ 

-pi 23 ) m - £ -^ 3 r^riVA f da a "" 1(1 - 
r(vi)r(v 2 )r(vz) j 



i 

db b m - £ -^- l (l - fym-e-^-l [!_„(!_ a .)]T»- e -^„ 





_/ 2 \m-e-i/i23 1 r{m — £ — V\)r(m — £ — 

~ ^ Pl2s) rWrfa) r(2m-2e- 1/123) 

~ r(n + ^)r(n-m + e + ^ a3 ) „ 

n— 

where x = VyiIVyizi ^23 = ^2 + ^3 and 1/123 = v \ + ^2 + ^3- To arrive at the 
last line of (J5J one expands [1 — a (1 — x)] m e_!yi23 according to 

r( C ) r tu r (- + 1 ) ' 
Then all Feynman parameter integrals are of the form 



1 





The infinite sum in the last line of (JSJ is a hypergeometric function, where 
the small parameter e occurs in the Gamma- functions. 

More complicated loop integrals yield additional classes of infinite sums. 
The following types of infinite sums occur: 

Type A: 

Y> rji + ax) r(i + a k ) i 

^r(i + ai)"T(i + a' fc ) 1 > 

Up to prefactors the hypergeometric functions j+\Fj fall into this class. 
The example discussed above is also contained in this class. 
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Type B: 

~ ~ r(j + ai ) r(i + ak ) r(j + h) r(j + b t ) 

hh r{i + a[) "' r{% + a ' k) r{j + b ' i] "' ni + b ' i) 

r(i+j + c\) r(i + j + c' m ) 

An example for a function of this type is given by the first Appell function 
F x . 
Type C: 



EE 

i=0 j=0 



i + r(i + ai) r(i + tifc) 



3 ) r(i + a[)"T(i + a' k ) 

r(i+j + ci) r(i + j + c 



r(i + j + c' 1 )~r(i + j + c' m ) 



x l y 3 (10) 



Here, an example is given by the Kampe de Feriet function S±. 
Type D: 



OO OO 

EE 

i=o j=o 



i + j\r(i + ai ) r(i + o fc ) r(j + h) r(j + b t ) 



j ) r(i + a [) rti + a'jrv + bti r(j + b\ 
r(i + j + Cl ) r(i+j + c. 



zy (ii) 



r(i+j + ci) r(i+j + d m 

An example for a function of this type is the second Appell function F% ■ 

Here, all a n , a' n , b ni b' n , c n and c' n are of the form "integer + const ■ e". The 
task is now to expand these functions systematically into a Laurent series in 
e, such that the resulting algorithms are suitable for an implementation into 
a symbolic computer algebra program. Due to the large size of intermediate 
expressions which occur in perturbative calculations, computer algebra sys- 
tems like FORM or GiNaC @|, which can handle large amounts of data 
are used. An implementation of the algorithms reviewed below can be found 
in0. 



3 Nested Sums 

In this section I review the underlying mathematical structure for the system- 
atic expansion of the functions in l|HJl-dJ- I discuss properties of particular 
forms of nested sums, which are called Z-sums and show that they form a 
Hopf algebra. This Hopf algebra admits as additional structures a conjugation 
and a convolution product. This summary is based on [B], where additional 
information can be found. Z-sums are defined by 
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Z(n;m 1 ,...,m k ;x 1 ,...,x k ) = V T^T • • • 7"^ ■ ( 12 ) 

n>ii>22>--->*fc>0 

fe is called the depth of the Z-sum and u> = mi + ... + is called the weight. 
If the sums go to Infinity (n = oo) the Z-sums are multiple polylogarithms 

Em 

Z(oo;m 1 , ...,m k ;x 1 , ...,x k ) = H mkt ... imi (x k , ...,x x ). (13) 

For x\ = ... = x k = 1 the definition reduces to the Euler-Zagier sums [7J|H]: 

Z(n;mi,...,TO fc ;l,...,l) = Z mij ... >mfc (n). (14) 

For 7i = oo and xi = ... = x& = 1 the sum is a multiple C- value 9 : 

Z(oo;mi, ...,m k ; 1, 1) = C(mk, ...,mi). (15) 

The multiple polylogarithms contain as the notation already suggests as sub- 
sets the classical polylogarithms Li„(x) lH, as well as Nielsen's generalized 
polylogarithms 

S n , P (x) = Li 1: ... !lin+1 (l, 1, x), (16) 

p-i 

and the harmonic polylogarithms |12| 

^mi,...,mfc(^) — Li mfc _ mi (l,...,l,x). (17) 

fc-i 

The usefulness of the Z-sums lies in the fact, that they interpolate between 
multiple polylogarithms and Euler-Zagier sums. 

In addition to Z-sums, it is sometimes useful to introduce as well S-sums. 
S-sums are defined by 

x h x ik 

S(n;mi,...,m k ;x 1 ,...,x k ) = 7^ • •• 7"^- ( 18 ) 

— *i *fe 

The S'-sums reduce for x x = ... = Xk = 1 (and positive rrij) to harmonic sums 
C3|: 

S(n;mi, ...,m k ; 1, 1) = S mi ,...,m k (n). (19) 

The S'-sums are closely related to the Z-sums, the difference being the upper 
summation boundary for the nested sums: (i — 1) for Z-sums, i for S-sums. 
The introduction of S-sums is redundant, since S-sums can be expressed in 
terms of Z-sums and vice versa. It is however convenient to introduce both 
Z-sums and S-sums, since some properties are more naturally expressed in 
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terms of Z-sums while others are more naturally expressed in terms of S'-sums. 
An algorithm for the conversion from Z-sums to S'-sums and vice versa can 
be found in 

The Z-sums form an algebra. The unit element in the algebra is given by 
the empty sum 

e = Z[n). (20) 

The empty sum Z(n) equals 1 for non-negative integer n. Before I discuss the 
multiplication rule, let me note that the basic building blocks of Z-sums are 
expressions of the form 

x" 

— > ( 21 ) 

which will be called "letters" . For fixed n, one can multiply two letters with 
the same n: 

e.g. the Xj's are multiplied and the degrees are added. Let us call the set of all 
letters the alphabet A. As a short-hand notation I will in the following denote 
a letter just by Xj = x T j /n mj . A word is an ordered sequence of letters, e.g. 

W = X 1 ,X 2 ,...,X k . (23) 

The word of length zero is denoted by e. The Z-sums defined in (|12|l are 
therefore completely specified by the upper summation limit n and a word 
W . A quasi-shuffle algebra A on the vectorspace of words is defined by [Tl] 

eoW = W oe = W, 
{X u W 1 ) o (X 2 , W 2 ) = X u (Wi o (X 2 , W 2 )) + X 2 , ((X u W x ) o W 2 ) 

+{X 1 -X 2 ),{W 1 oW 2 ). (24) 

Note that "•" denotes multiplication of letters as defined in eq. (|2*2|) . whereas 
"o" denotes the product in the algebra A, recursively defined in eq. (|2*4*|l . This 
defines a quasi-shuffle product for Z-sums. The recursive definition in 124|) 
translates for Z-sums into 

Zmi^.^rrik ij^) * Z m / _ _ _ [jl) 
n 1 

= ^ y JmT^m2,...,mk (*1 — ^-)Z m ' 1 ,...,m' l {}l ~ 1) 



I 

i 1= l 



1 

1 



■ X! ~^T Z ™u-;m k ('i'2 - l)Z m ' 2 ,...,m' l (i2 - 1) 
i 2 = l l 2 



n 1 

ZZ jmi+mi Z m»,...,mS ~ 1 )Z m '^...,m' l {i ~ !)• (25) 



i=! 1 
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1 2 



• • • • 



• • • • 



1-1 



1 2 



1 1 



h 



h 



Fig. 2. Sketch of the proof for the multiplication of Z-sums. The sum over the 
square is replaced by the sum over the three regions on the r.h.s. 



The proof that Z-sums obey the quasi-shuffle algebra is sketched in Fig. [21 
The outermost sums of the Z-sums on the l.h.s of (|25f) are split into the three 
regions indicated in Fig. |3 A simple example for the multiplication of two 
Z-sums is 



Z(n;mi;xi)Z(n;mi;x2) = 
Z(n; mi, m 2 ; Xi, x^) + Z(n; m 2 , m±; x 2l xi) + Z(n; m\ + m 2 ; xix 2 ). (26) 

The quasi-shuffle algebra A is isomorphic to the free polynomial algebra 
on the Lyndon words. If one introduces a lexicographic ordering on the letters 
of the alphabet A, a Lyndon word is defined by the property 



W < V 



(27) 



for any subwords U and V such that W = U, V. Here U, V means just con- 
catenation of U and V. 

The Z-sums form actuall a Hopf algebra. It is convenient to phrase the 
coalgebra structure in terms of rooted trees. Z-sums can be represented as 
rooted trees without any sidebranchings. As a concrete example the pictorial 
representation of a sum of depth three reads: 



Z(n; mi, m 2 , m 3 ; xi,x 2 ,x 3 ) 



EEE — 

»i=l »a=l 13=1 



H mi i 2 m2 *3 m3 



X\ f 

X2 
X3 * 



(28) 



The outermost sum corresponds to the root. By convention, the root is al- 
ways drawn on the top. Trees with sidebranchings are given by nested sums 
with more than one subsum, for example: 



E ^r z ^ ~ 1; m2 ' x ^) z ( i - !; m s; x 3) 
»=i 



■r 2 




(29) 



X 3 - 
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Of course, due to the multiplication formula, trees with sidebranchings can 
always be reduced to trees without any sidebranchings. The coalgebra struc- 
ture is now formulated in terms of rooted trees. I first introduce some notation 
how to manipulate rooted trees, following the notation of Kreimer and Connes 
An elementary cut of a rooted tree is a cut at a single chosen edge. An 
admissible cut is any assignment of elementary cuts to a rooted tree such that 
any path from any vertex of the tree to the root has at most one elementary 
cut. An admissible cut maps a tree t to a monomial in trees t\ o ... o tk+i- 
Note that precisely one of these subtrees tj will contain the root of t. Denote 
this distinguished tree by R c (t), and the monomial delivered by the k other 
factors by P (t). The counit e is given by 

e(e) = 1, 

g(t) = 0, t ^ e. (30) 
The coproduct A is defined by the equations 
A(e) = e <8> e, 

A(t) = e®t + t®e+ P C (t)<Z)R c (t), 

adm. cuts O of t 

A(h o ... o t k ) = Z\(ii)(o ® o)...(o ® o)A(t k ). (31) 
The antipode S is given by 

S(e) = e, 

S(t)=~t- E S(P c (t))oR c (t), 

adm. cuts C Of t 

S(h o ... o t k ) = S{tx) o ... o S(t k ). (32) 

Since the multiplication in the algebra is commutative the antipode satisfies 

S 2 = id. (33) 

Let me give some examples for the coproduct and the antipode for Z-sums: 

AZ(n; mi; x%) = e® Z(n; mi; Xi) + Z(n; mi; Xi) ® e, 
AZ(n; mi, ma; 2:1,12) = e <8> i?(n; mi, m2; 2:1, 2:2) + Z(ra; mi, m2; 2:1,2:2) ® e 
+Z(n; m 2 ; x 2 ) ® mi; xi), (34) 

SZ(n; mi; 2:1) = —Z(n;mi;xi), 
SZ(n; mi,m 2 ; 2:1,2:2) = Z(n; m 2 , mr,X2, xi) + Z(n; mi + mr, 2:12:2). (35) 

The Hopf algebra of nested sums has additional structures if we allow expres- 
sions of the form 
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'o 



n" 



Z(n\m Xl ...,m k ;xi, ...,x k ), 



(36) 



e.g. Z-sums multiplied by a letter. Then the following convolution product 

n—l ,; 



y =-z(i-i\ 



v 



(n-i) 



- — 7 Z(n — i — 1; ...) 



(37) 



can again be expressed in terms of expressions of the form H3(jf) . An example 



is 



E -zi^-i) 
^ — ^ i 



n 



i=l 

a; 
n 



(n-i) 



—Z\(n — i — 1) 



v i i i 1 i y x y 

Z I n- 1; 1,1,1; -, -, - 

x y x 



Z\ n-1; 1,1,1;-, 1, - 

a; y 



Z(n-1; 1,1,1; 1,^,1 



+ <-> y) . 



In addition there is a conjugation, e.g. sums of the form 



E 



(38) 



(39) 



can also be reduced to terms of the form l|3tj|) . Although one can easily con- 
vert between the notations for .S'-sums and Z-sums, expressions involving a 
conjugation tend to be shorter when expressed in terms of S'-sums. The name 
conjugation stems from the following fact: To any function f(n) of an integer 
variable n one can define a conjugated function Co f(n) as the following sum 



Co/(n) = E(" ) ("!)'/( 



(40) 



Then conjugation satisfies the following two properties: 

Col = 1, 
CoCof(n)=f(n). 

An example for a sum involving a conjugation is 



(41) 



E " Wt^o 



5* ( n; 1, 1; 1 — x, 



1 - x 



S" (n; 1, 1; 1 — a;, 1) . 



(42) 



Finally there is the combination of conjugation and convolution, e.g. sums of 
the form 
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£ ? (-1)' 



(n — i) 1 



-t-S^tj - i; ...) 



(43) 



can also be reduced to terms of the form . An example is given by 



-E 1 (^)(- 1 ) 1 S(i;l;x)S(n-i;l;y) = 

i=l \ ' 



-<^(n;l;y) + (l-z) 7 
n 



S [n; 1; 



1 



- 5 n: 1 



(-1)' 



S(n;l;x) + (!-</)" 



5 n; 1; 



S\n;l\- \ 



(44) 



4 Expansion of hypergeometric functions 

In this section I discuss how the algebraic tools introduced in the previous 
section can be used to solve the problems outlined at the end of Sect. [3 First 
I give some motivation for the introduction of Z-sums: The essential point is 
that Z-sums interpolate between multiple polylogarithms and Euler-Zagier- 
sums, such that the interpolation is compatible with the algebra structure. 
On the one hand, we expect multiple polylogarithm to appear in the Laurent 
expansion of the transcendental functions Q- ljTJ , a fact which is confirmed 
a posteriori. Therefore it is important that multiple polylogarithms are con- 
tained in the class of Z-sums. On the other the expansion parameter e occurs 
in the functions l!51)- (lll|) inside the arguments of Gamma- functions. The basic 
formula for the expansion of Gamma-functions reads 

r(n + e)= T(l + e)r{n) [l + eZ x {n - 1) + e 2 Z n {n - 1) 

+e 3 Z nl (n - 1) + ... + e n - 1 Z n ... 1 (n - 1)] , (45) 

containing Euler-Zagier sums for finite n. As a simple example I discuss the 
expansion of 

E r(i + ai + tis)r(i + a 2 + t 2 e) i , . 

n r(i + i)r(i + a 3 + t 3 e) x [ ' 

into a Laurent series in e. Here oi, a 2 and 03 are assumed to be integers. Up 
to prefactors the expression in l|4t)|l is a hypergeometric function 2 Fi. Using 
r(x + 1) = xT(x), partial fractioning and an adjustment of the summation 
index one can transform 146(1 into terms of the form 

^ r(i)r(i + t 3 s) 1 ' 
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where m is an integer. Now using l|45|) one obtains 

r(1 + £) |t (1 + e^ (*-!) + ...) ^ (48) 

Inverting the power series in the denominator and truncating in e one obtains 
in each order in e terms of the form 



^ —Z mi ... mk (i - tyZm'^.mfo - l)Z m //... m »(z - 1) (49) 



i=l 



Using the quasi-shufne product for Z-sums the three Euler-Zagier sums can 
be reduced to single Euler-Zagier sums and one finally arrives at terms of the 
form 

oo i 

^2^Z mi ... mk (i-l), (50) 

i=l 

which are harmonic polylogarithms H mmii „ lk (x). This completes the al- 
gorithm for the expansion in e for sums of the form (|SJ|. Since the one- loop 
integral discussed in JSJ is a special case of ||SJ|, this algorithm also applies to 
the integral JSJ. In addition, this algorithm shows that in the expansion of 
hypergeometric functions j+iFj(ai, aj+i; bi, bj; x) around integer val- 
ues of the parameters a>k and bi only harmonic polylogarithms appear in the 
result. 

The algorithm for the expansion of sums of type QB} used the multiplication 
formula for Z-sums to pass from (|49|) to l|50[) . To expand double sums of type 
(0 one needs in addition the convolution product To expand sums of 

type qi0[) the conjugation l|39fl is needed. Finally, for sums of type i|ll|) the 
combination of conjugation and convolution as in (|43|l is required. More details 
can be found in [5]. 

Let me come back to the example of the one-loop Feynman integral dis- 
cussed in Sect. EJ For v\ = v-i = = 1 and m — 2 in (0 one obtains: 



/ 



1 1 1 



r(-e)r(i e)r(i + £ ) (-p^)" 1 " 6 f> ^ „ n , 

= W^M) ^~tx ^ 

Here, all harmonic polylogarithms can be expressed in terms of Nielsen poly- 
logarithms, which in turn simplify to powers of the standard logarithm: 

■x) =S , n (l-x) = ^-(]nx) n . (52) 
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This particular example is very simple and one recovers the well-known all- 
order result 



r(i-e) 2 r(i + e) (-p? 23 ) 1 E 1 - x 



(53) 



r(l-2e) e 2 l-x ' 

which (for this simple example) can also be obtained by direct integration. 

5 Multiple polylogarithms 

The multiple polylogarithms are special cases of Z-sums. They are obtained 
from Z-sums by taking the outermost sum to infinity: 

Z(ao;mi, —,m k ;xx, ...,x k ) = hi mki ... iTni (x k , ...,x%). (54) 

The reversed order of the arguments and indices on the r.h.s. follows the nota- 
tion of Goncharov |23j . They have been studied extensively in the literature by 
physicists |12l I13| . J7j-|2J| and mathematicians [§]>|22]~|I12|- Here I summa- 
rize the most important properties. Being special cases of Z-sums they obey 
the quasi-shuffle Hopf algebra for Z-sums. Multiple polylogarithms have been 
defined in this article via the sum representation (|T3|) . In addition, they admit 
an integral representation. From this integral representation a second algebra 
structure arises, which turns out to be a shuffle Hopf algebra. To discuss this 
second Hopf algebra it is convenient to introduce for z k ^ the following 
functions 

G(*x,..,*;v)= fr^- / * ./"-J**-. (55) 
J h - zi J t 2 - z 2 J tk - z k 



In this definition one variable is redundant due to the following scaling rela- 
tion: 

G(zi,...,z k ;y) = G{xz\, ...,xz k ;xy) (56) 

If one further defines 

g(z;y) = —^— , (57) 

y-z 



then one has 



and 



^G{zi,...,z k ;y) = g(z 1 ;y)G{z2,...,z k ;y) (58) 
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z 1 ,z 2 ,...,z k ;y) = J dt g(z 1 ;i)G(z 2 , z k ;t). (59) 



G( 



One can sligthly enlarge the set and define G(0, 0; y) with k zeros for z\ to 
z k to be 

G(0,...,0;y) = l(lny) fe . (60) 

This permits us to allow trailing zeros in the sequence (z±, z k ) by defining 
the function G with trailing zeros via Ij59(l and H6()|l . To relate the multiple 
polylogarithms to the functions G it is convenient to introduce the following 
short-hand notation: 

Gm u ...,m h (zi, —,Zk;y) = G(O^0, Z\ , .. ., z k -i , 0^0, z k ; y) (61) 

mi— 1 m k ~l 

Here, all Zj for j = 1, fc are assumed to be non-zero. One then finds 

Lijrfcfc,...,mi { x k> ■■■■> x l) — (~l-) k G mi ,...,m k ( — j ) - --j i 1 I ■ (62) 



X\ X\Xi X\...X k 



The inverse formula reads 



G mi ,...,m k {zi, -,z k ;y) = (-l) fc Li mfe> ... >mi ( — , — ) . (63) 

V ^fe z 2 ZiJ 

Eq. Ij62(l together with (|61|l and (|55(l defines an integral representation for the 
multiple polylogarithms. To make this more explicit I first introduce some 
notation for iterated integrals 

A A t n t 2 

dt dt f dt n f dt n - X f dt x 

o ... o = / / x ... x / (64) 

t — a n t — a\ J t n — a n J — a„_i J t\ — a% 

o oo o 

and the short hand notation: 

A A 

dt \ m dt f dt dt dt 

— o = 0...— 0- . (65) 

t J t — a J t t t — a 
o v ^ ' 

m times 

The integral representation for 'Li m , k ,...,m 1 (x k , ...,xi) reads then 

f ^ < ukf (dt \ mi -' dt 

Li mjfc ,...,mi(a;fe, —,xi) = (-1) / — o 



t J t-bi 
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where the bj 's are related to the Xj 's 

bj = . (67) 

X\X%...Xj 

From the iterated integral representation (|55l) a second algebra structure for 
the functions G(zi, z k ;y) (and through (|6~2)l also for the multiple polyloga- 
rithms) is obtained as follows: We take the Zj's as letters and call a sequence 
of ordered letters w — Zi,...,z k a word. Then the function G(zi, ...,Zf.',y) is 
uniquely specified by the word w = z\, z k and the variable y. The neutral 
element e is given by the empty word, equivalent to 

G(;y) = l. (68) 

A shuffle algebra on the vector space of words is defined by 

eow = woe = w, 
{zi,w\) o (z 2 ,w 2 ) = zi, (wi o (z 2 , w 2 )) + z 2 , ((zi, wi) o w 2 ). (69) 

Note that this definition is very similar to the definition of the quasi-shufne 
algebra H24fl . except that the third term in (|24|l is missing. In fact, a shuffle 
algebra is a special case of a quasi-shuffle algebra, where the product of two 
letters is degenerate: Xi ■ X 2 = for all letters X\ and X 2 in the notation of 
Sect.|3 The definition of the shuffle product Ij69(l translates into the following 
recursive definition of the product of two G- functions: 

G(zi, z k ;y) x G(z k+1 , z n ; y) = (70) 
v 

G(z 2 , z k ;t)G(z k+ i, z n ; t) 

t — Zl 

o 

y 

G(zi, ...,z k ;t)G(z k+2l ...,z n ;t) (71) 

t — z k+ i 
o 

For the discussion of the coalgebra part for the functions G(zi, ...,z k ;y) 
we may proceed as in Sect. |21 and associate to any function G(zi, ...,z k ;y) a 
rooted tree without sidebranchings as in the following example: 



Zl 

G(zi,Z2,z 3 ;y) = z 2 

Z3 



(72) 



The outermost integration (involving zi) corresponds to the root. The for- 
mulae for the coproduct (|31() and the antipode l|32l) apply then also to the 
functions G(zi, ...,z k ;y). 



16 Stefan Weinzierl 

A shuffle algebra is simpler than a quasi-shufffe algebra and one finds for a 
shuffle algebra besides the recursive definitions of the product, the coproduct 
and the antipode also closed formulae for these operations. For the product 
one has 

G(z 1 ,...,z k ;y) G(z k+1 ,...,z k+ i;y) = ^2 G ( z <x(i)' -> z cr(k+i)', y) , (73) 

shuf fie 

where the sum is over all permutations which preserve the relative order of the 
strings z%, z k and z k +\, z k +i- This explains the name "shuffle product" . 
For the coproduct one has 

k 

AG(zi,...,z k ;y) = G(zi, Zf, y) ® G(z j+1 , z k ; y) (74) 

3=0 

and for the antipode one finds 

SG( Zl , z k -y) = (-l) k G(z k , Zl ;y). (75) 
The shuffle multiplication is commutative and the antipode satisfies therefore 

S 2 = id. (76) 

From H75f) this is evident. 



6 The antipode and integration-by-parts 

Integration-by-parts has always been a powerful tool for calculations in parti- 
cle physics. By using integration-by-parts one may obtain an identity between 
various G-functions. The starting point is as follows: 

v 

d 



G(zi,...,z k ;y) = / dt ( — G{z x ;t) ) G(z 2 ,...,z k ;y) 



o 

y 



G(z 1 ;y)G(z 2 ,...,z k ;y) - / dt G(z 1 ;t)g(z 2 ;t)G(z 3 , z k ;y) 



o 

y 

= G(z 1 ;y)G(z 2 ,...,z k ;y)- J dt (J^G(z 2 , z^t^J G(z 3 , z k ;y). (77) 

o 

Repeating this procedure one arrives at the following integration-by-parts 
identity: 

G{zi,...,z k ;y) + {-l) k G(z kl Zl ; y) 
= G(z 1 ;y)G(z 2 , ...,z k ;y) - G(z 2 , z 1 ;y)G(z 3 , ...,z k ;y) + ... 

-{-l) k - 1 G(z k - 1 ,...zi;y)G(z k ;y), (78) 
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which relates the combination G(z\, z k ; y) + (—l) k G(z k ,...,zi;y) to G- 
functions of lower depth. This relation is useful in simplifying expressions. 
Eq. H78|) can also be derived in a different way. In a Hopf algebra we have for 
any non-trivial element w the following relation involving the antipode: 

-S{wW) = 0. (79) 

O) 

Here Sweedler's notation has been used. Sweedler's notation writes the co- 
product of an element w as 

A ( w ) = z2 w{1) ® w(2) - ( 80 ) 

(w) 

Working out the relation l|79|) for the shuffle algebra of the functions G(zx, 
Zk',y), we recover (|7S|) . 

We may now proceed and check if (|79|l provides also a non-trivial relation 
for the quasi-shufne algebra of Z-sums. This requires first some notation: A 
composition of a positive integer k is a sequence / = (ix,...,ii) of positive 
integers such that i\ + ...ii = k. The set of all composition of k is denoted by 
C(k). Compositions act on Z-sums as 

(ii, o 2(fi;mi, ...,m fe ;xi, ...,x k ) 

= Z (n; mi + ... + m ii; m il+ x + ... + rn il+i2 , m il+ ... +il _ 1+ i + ... 

+ j X\...Xi^ , Xi 1 J r \...Xi 1 -\-i 2 , Xi 1 J r ,,,-\^i l l J r i...Xi 1 ^_^i l ) , (81) 

e.g. the first i\ letters of the Z-sma are combined into one new letter, the next 
i2 letters are combined into the second new letter, etc.. With this notation for 
compositions one obtains the following closed formula for the antipode in the 
quasi-shufne algebra: 

SZ(n;mi, ...,m k ;xi, ...,x k ) = (-l) fe 1 ° z ( n : m k, ■-,rn 1 ;x k , ...,xi) 

iec(k) 

(82) 

From 179|l we then obtain 

Z(n;mx, ...,m k ;xx, ...,x k ) + (-l) k Z(n; m k , ...,mi;x k , -,Xx) 
= - P C (Z(n;mi,...,m k ;x 1 ,... 1 x k )) 

adm. cuts 

■S (R c (Z(n;mi, ...,m k ;xi, ...,x k ))) 

-(- l ) k X! Io Z(n;m k ,...,mi;x k ,...,xi). (83) 
zec(fe)\(i,i,...,i) 

Again, the combination Z(n; mi, m k ; X\, x k ) + (—l) k Z(n; m k , mi; x k , 
...,xi) reduces to Z-sums of lower depth, similar to (|78|l . We therefore ob- 
tained an "integration-by-parts" identity for objects, which don't have an 
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integral representation. We first observed, that for the G-functions, which 
have an integral representation, the integration-by-parts identites are equal 
to the identities obtained from the antipode. After this abstraction towards 
an algebraic formulation, one can translate these relations to cases, which only 
have the appropriate algebra structure, but not necessarily a concrete integral 
representation. As an example we have 

Z(n; mi, m 2 , m 3 ; xi, x 2 , x 3 ) - Z{n; m 3 , m 2 , mi; S3, x 2 , xi) = 
Z(n\ m\\x\)Z(n; m 2 , m 3 ; x 2 , x 3 ) - Z(n; m 2 , mi\x 2 ,xi)Z{n\ m 3 ; x 3 ) 
—Z(n; mi + m 2 ; xix 2 )Z(n; m 3 ; x 3 ) + Z(n; m 2 + m 3 , mi; x 2 x 3 , X\) 
+Z(n; m 3 , mi + m 2 ; x 3 , X\x 2 ) + Z(n; mi + m 2 + m 3 ;xix 2 x 3 ), (84) 

which expresses the combination of the two Z-sums of depth 3 as Z-sums 
of lower depth. The analog example for the shuffle algebra of the G-function 
reads: 

G(zi,z 2 ,z 3 ;y) - G(z 3 ,z 2 ,z 1 ;y) = G(z 1 ;y)G(z 2 ,z 3 ;y) - G(z 2 ,z 1 ;y)G(z 3 ;y). 

(85) 

Multiple polylogarithms obey both the quasi-shufhe algebra and the shuffle 
algebra. Therefore we have for multiple polylogarithms two relations, which 
are in general independent. 

7 Summary 

In this article I discussed the mathematics underlying the calculation of Feyn- 
man loop integrals. The algorithms are based on Z-sums, which form a Hopf 
algebra with a quasi-shufHe product. This algebra has as additional structure 
a conjugation and a convolution product. In the final results multiple poly- 
logarithms appear. Multiple polylogarithms obey apart from the quasi-shufHe 
algebra a second Hopf algebra. This additional Hopf algebra has a shuffle 
product. 

A Notations and conventions 

There are several notations for the multiple polylogarithms. I briefly summ- 
marize them here. In this article multiple polylogarithms arc defined via the 
sum representation 

Lim fc ,...,mi(^fc, = Z(oo;mi, ■■■,m k ;xi, ...,x k ). (86) 

The reversed order of the arguments and indices for Li mk> .., imi (xk, --tXi) fol- 
lows the notation of Goncharov [231 ■ Gehrmann and Remiddi |181 1191 12U] use 
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the notation G(z\, ...,Zk',y) and G mi) ... jjnji (z[, z' k ; y). The relation with the 
notation above is 

Li mk ,..., mi (xk, = (-l) fc Gm 1 ,...,m fe (^1,&2, -, &fc5 1) , (87) 

where 

b, = . (88) 

Borwein, Bradley, Broadhurst and Lisonek 9 denote multiple polylogarithms 

as 

/ uii,...,mk\ j • , \ , Qrn 

b bk J = ^hn k ,..., mi {Xk,-,Xl) (89) 

In the french literature HOI harmonic polylogarithms are often denoted as 

,...,mk \X) — ^mx,---,fnk\X) (90) 

and referred to as "multiple polylogarithms of a single variable". Note the 
order of the indices for Li mii ... iTrHi (x) in (|9U|) . 



References 

1. O. V. Tarasov, Phys. Rev. D54, 6479 (1996), hep-th/9606018 

2. O. V. Tarasov, Nucl. Phys. B502, 455 (1997), hep-ph/9703319 

3. J. A. M. Vermaseren, [math-ph/0010025| 

4. C. Bauer, A. Frink, and R. Kreckel, J. Symbolic Computation 33, 1 (2002), 
Ics.se/00040i5l 

5. S. Weinzierl, Comput. Phys. Commun. 145, 357 (2002), math-ph/0201011 

6. S. Moch, P. Uw er, and S. Weinzierl, J. Math. Phys. 43, 3363 (2002), 
|hep-ph/01 10083] 

7. L. Euler, Novi Comm. Acad. Sci. Petropol. 20, 140 (1775). 

8. D. Zagier, First European Congress of Mathematics, Vol. II, Birkhauser, Boston, 
497 (1994). 

9. J. M. Borwein, D. M. Bradley, D. J. Broadhurst and P. Lisonek, Trans. Amer. 
Math. Soc. 353:3, 907 (2001), math.CA/9910045 

10. L. Lewin, "Polylogarithms and associated functions", (North Holland, Amster- 
dam, 1981). 

11. N. Nielsen, Nova Acta Leopoldina (Halle) 90, 123 (1909). 

12. E. Remiddi and J. A. M. Vermaseren, Int. J. Mod. Phys. A15. 725 (2000). 
hep-ph/9905237 

13. J. A. M. Vermaseren, Int. J. Mod. Phys. A14, 2037 (1999), hep-ph/9806280 

14. M. E. Hoffman, J. Algebraic Combin. 11, 49 (2000), math.QA/9907173 

15. D. Kreimer, Adv. Theor. Math. Phys. 2, 303 (1998), q-alg/9707029 

16. A. Connes and D. Kreimer. Commun. Math. Phys. 199. 203 (1998). 
|hep-th/9808042| 



20 Stefan Weinzierl 



17. T. Gehrmann and E. Remiddi, Nucl. Phys. B601, 248 (2001), hep-ph/0008287 

18. T. Gehrmann an d E. Remiddi, Comput. Phys. Commun. 141, 296 (2001), 
|hep-ph/0107173| 

19. T. Gehrmann an d E. Remiddi, Comput. Phys. Commun. 144, 200 (2002), 
|hep-ph/01lT255] 

20. T. Gehrmann and E. Remiddi, Nucl. Phys. B640, 379 (2002), hep-ph/0207020 

21. S. Moch, P. U wer, and S. Weinzierl, Phys. Rev. D66, 114001 (2002), 
|hep-ph/0207043l 

22. R. M. Hain, alg-geom/9202022 

23. A. B. Goncharov, Math. Res. Lett. 5 , 497 (1998). 

24. A. B. Goncharov, | math.AG/0103059 | 

25. A. B. Goncharov, math. A G/0207036 

26. A. B. Goncharov, math.AG/0208144 

27. P. Elbaz-Vincent and H. Gangl, Comp. Math. 130, 161 (2002), 
|math.KT/0008089| 

28. H. Gangl, |math .KT/0207222 

29. H. M. Minh, M. Petitot and J. van der Hoeven, Discrete Math. 225:1-3, 217 
(2000). 

30. P. Cartier, Seminaire Bourbaki , 885 (2001), in french. 

31. J. Ecalle, Preprint Orsay 2002-23, (2002), in french with additional grammatical 
inventions, http:/ /www. math. u-psud.fr/~biblio/ppo/2002/ppo2002-23. html 

32. G. Racinetr |math.QA/0202142| in french. 



